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ABSTRACT 

A new model for cooling flows in X-ray clusters, capable of naturally 
explaining salient features observed, is proposed. The only requirement is that 
a significant relativistic component, in the form of cosmic rays (CR), be present 
in the intra-cluster medium and significantly frozen to the thermal gas. Such 
an addition qualitatively alters the conventional isobaric thermal instability 
criterion such that a fluid parcel becomes thermally stable when its thermal 
pressure drops below a threshold fraction of its CR pressure. Consequently, the 
lowest possible temperature at any radius is about one third of the ambient 
temperature at that radius, exactly as observed, In addition, we suggest that 
dissipation of internal gravity waves, excited by radial oscillatory motions of 
inward drifting cooling clouds about their radial equilibrium positions, may 
be responsible for heating up cooling gas. With the ultimate energy source 
for powering the cooling X-ray luminosity and heating up cooling gas being 
gravitational due to inward drifting cooling clouds as well as the general inward 
flow, heating is spatially distributed and energetically matched with cooling. 
One desirable property of this heating mechanism is that heating energy is 
strongly centrally concentrated, providing the required heating for emission-line 
nebulae. 



Subject headings: Clusters of galaxies - cooling flows - cosmic rays 



1. Introduction 

It is a firmly established observational fact that a substantial central region of many 
X-ray clusters of galaxies is cooling in X-rays on time scales that are shorter than their 
presumed age, implying that a large amount of gas would cool to lower temperatures, if not 
intervened (e.g., Fabian 1994; White, Jones, & Forman 1997; Peres et al. 1998; Allen 2000). 
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The classic picture was that the gas would cool isobarically and sink to the center 
(Fabian & Nulsen 1977; Cowie & Binney 1977; Fabian 1994). While simple and elegant, 
this model is at odds with current observational evidence, with two major contradictions. 
First, the amount of cooled gas at the centers of clusters, although having been detected 
in various forms, including UV emitting gas (e.g., Oegerle et al. 2001), Ha emitting gas 
(e.g., Heckman et al. 1989; Crawford et al. 1999), molecular gas (e.g., Jaffe & Bremer 1997; 
O'Dea et al. 1998; Donahue et al. 2000 Edge 2001; Salome & Combes 2003; Edge & Frayer 
2003) and star formation (e.g., Fabian et al. 1991; Hansen, Jorgensen, & Norgaard-Nielsen 
1995; Allen 1995; Smith et al. 1997), appears to be far short of what is implied by the 
cooling rate in X-rays. This is the "classic cooling- flow problem" (Peterson et al. 2003). 
Second, the supposedly cooled X-ray gas does not appear to show up at temperatures lower 
than about one third of the maximum temperature at each radius (e.g., Peterson et al. 
2001,2003; Kaastra et al. 2004); i.e., the expected soft X-rays are missing (e.g., Fabian et 
al. 2002) — the "soft X-ray cooling-flow problem" (Peterson et al. 2003). 

In light of these observational developments various amendments to the standard 
cooling flow model with isobaric thermal instability have been proposed to change the state 
of the X-ray cooled gas based on energetic considerations. There are mostly two schools 
of ideas, either with heating by some sources or by cooling in channels other than soft 
X-rays. The reader is referred to Mathews & Brighenti (2003) and Peterson et al. (2003) 
for up-to-date assessments of various proposed models, confronted with observations. In 
short, it appears that none of the proposed amendments is fully satisfactory. 

In this Letter we put forth a new model, based on the assumption that a significant 
fraction of total pressure of the X-ray cluster gas is in the form of cosmic rays, an 
unavoidable by-product formed at the ubiquitous high velocity shocks during the formation 
of X-ray clusters. We will show that this model explains all major characteristics of 
observed X-ray cooling flows naturally. 

2. A Two-Fluid Model 

Our purpose is to present a viable physical model for the cooling flows in X-ray clusters. 
Rigorous calculations possibly require detailed simulations and are not attempted here. 
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2.1. Cosmic-Ray Pressure in the Intra-Cluster Medium 



Astrophysical shocks have long been recognized as efficient sites for CR acceleration 
(Blandford & Ostriker 1978). Clusters of galaxies are full of shocks under the action of 
gravity, including accretion shocks (e.g., Bertschinger 1985) and internal shocks (e.g., 
Miniati et al. 2000), and hence should be excellent production sites of CR. Miniati et al. 
(2001), through detailed shock acceleration process of CR, find that shocks, especially 
internal high energy (due to high density and high velocity) shocks within clusters can 
convert a large fraction of gravitational energy released during the formation of X-ray 
clusters into CR. Miniati et al. (2001) show that a few tens of percent of the total pressure 
in the intra-cluster medium (ICM) may be in the form of CR, while Ryu et al. (2003) 
subsequently show that about 1/3 of the total pressure in the ICM may be in CR. We 
define a = Pcii/Pt, where Pcr and P t are the CR and total pressure, respectively, and for 
simplicity a is assumed to be constant spatially. It is assumed that thermal gas is coupled 
on the scales concerned to CR by magnetic field. Significant magnetic field is known to 
exist in clusters of galaxies (e.g., Kim et al. 1990; Taylor & Perley 1993; Feretti et al. 1995), 
although for the purpose of coupling CR to thermal gas a strong magnetic field is not 
required. Observationally, evidence for the existence of a significant CR pressure in clusters 
of galaxies is also strong (e.g., Blasi 1999). In passing, we note that our conclusions will 
be the same, if CR pressure is replaced by another form of relativistic component such as 
tangled magnetic fields. 



Consider a somewhat overdense parcel of gas containing a significant fraction of CR 
pressure in X-ray clusters. Such a parcel would normally tend to fall towards the center, 
because of anti-buoyancy and thermal instability. We show that, if CR are frozen to 
the thermal gas in the parcel, the thermal gas becomes thermally stable below a certain 
temperature. Assuming the ratio of specific heats for CR to be 4/3 and denoting the ratio 
of the thermal gas pressure to the CR pressure in the parcel as (3, then, the standard (local) 
isobaric thermal instability criterion (Field 1965; Balbus 1985) can be shown to become 



2.2. Dynamics of Cosmic Ray Modulated Cooling Flows 




(1) 



Here ( = ApT 1 / 2 — r](p) is the net cooling rate per unit mass, where the first term is the 
Bremsstrahlung cooling with A being a constant and the second term rj is the total heating 
term, which may be a function of many things. We will consider two simple cases: 1) i] is 
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constant, meaning that heating is uniform per unit mass; 2) r\ oc p™ 1 , meaning that heating 
is uniform per unit volume. It seems plausible that these two cases may bracket possible 
real situations. The following analysis may be straight-forwardly extended to more complex 
heating functions which may also depend on gas temperature. Inserting ( for the first case 
with equal mass heating into equation (1), the thermal instability criterion simplifies to 

P > |. P) 

For the second case with uniform volume heating we make the following assumption in 
order to allow further progress: there is no net cooling at the ambient gas at each radius 
to normalize the heating term, i.e., ( = ApT 1 ! 2 — Ap 2 T^ 2 p -1 , where p a and T a are the 
ambient gas density and temperature, respectively; with this form equation (1) becomes 

4 

P > 3(l + 2(£)-'({)-V2)- (3) 

In both cases we see that the conventional picture of thermal instability in X-ray cluster 
cooling flows (Mathews & Bregman 1978; Fabian, Nulsen, & Canizares 1984; Nulsen 1986) 
applies only in a limited sense, i.e., within a finite temperature range. But, in the presence 
of CR pressure, ultimately, there will be no runaway thermal instability and a temperature 
floor is put in place and further cooling of clouds relative to the ambient gas ceases. 

In the finite range of temperature where thermal instability does apply, the denser gas 
parcel cools, gets compressed and tends to fall inward. Since CR virtually do not cool, 
the CR pressure increases relatively during this infall process. It would finally reach a 
quasi-equilibrium position at a smaller radius in the comoving flow, with the exact location 
depending on the density and temperature run in the cooling flow region. Two conditions 
are met there. First, the gas parcel is in dynamical equilibrium at that radius which 
implies p p = p a , where p p and p a are the gas density in the parcel and ambient gas density, 
respectively. Second, the parcel needs to be in pressure equilibrium with the ambient gas, 
requiring p p = p a , where p p — (1 + 1/ j3)p p kT p /m p and p a = j^p a kT a /m p are the total 
pressure in the parcel and in the ambient gas, respectively, k is the Boltzmann's constant, 
m p is the proton mass. Combining these relations gives 

Tp = (1 - a )(l + P) Ta ' (4) 

Thus, a floor value for (3 for thermal instability translates to a floor in temperature. We 
find that for the case with constant heating per unit mass the lowest possible temperature 
for a cooling gas parcel is 

4 1 
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whereas for the case with constant heating per unit volume (and utilizing dynamical 
equilibrium condition p = p a ), we obtain 



18 + ^-6,/9 + ^ 

rp 1 — a V 1—a rp lp\ 

I P~ 4g J-a- \0) 

Miniati et al. (2001) show that a is about a few tens of percent, while Ryu et al. (2003) show 
a ~ 1/3. Thus, using their predicted range a = 10 — 30% we obtain the floor temperature 

T p = (0.6 - 0.8)T a , (7) 

if heating is uniform in mass, and 

T p = (0.2 - 0.3)T a , (8) 

if heating is uniform in volume. This indicates that the temperature floor may somewhat 
vary, depending on detailed dependence of heating on density. But the essential point is 
that there will always be an absolute floor at about one third of the ambient temperature 
at any radius. Thus, it is borne out in our model, without any fine tuning, that the coldest 
gas at each radius that coexists with the hotter ambient gas may have a temperature of 
about one third of that of the ambient gas at that radius, explaining the seemingly puzzling 
observations (Peterson et al. 2003; Kaastra et al. 2004). 

Equation (8) only demands that the ratio of the coldest temperature to ambient 
temperature at any radius maintain a certain ratio. However, it does not prevent overall 
ambient temperature T a from decreasing or increasing. In the cooling flows T a would 
decrease in the absence of a balanced heating source. In the next section we suggest a 
heating mechanism that is physically based and has the desired properties. 



2.3. Heating by Dissipation of Internal Gravity Waves 

The rising gas density and dropping gas temperature towards the center of cooling flows 
produce a subadiabatic run and therefore gas is convectively stable under the Schwarschild 
criterion for adiabatic gas. However, for gas parcels with temperature ranging from ~ l/3T a 
to T a , we have a thermally unstable gas, which is also convectively unstable even for a 
subadiabatic run (Defouw 1970). But for a strongly subadiabatic run, as in the cooling 
flow regions where cooling time may be longer than the inverse of oscillation frequency, 
it is expected that the thermal-convective instability actually becomes overstable, first 
suggested analytically (Defouw 1970) and subsequently confirmed by numerical simulations 
(Malagoli, Rosner, & Fryxell 1990), manifested by the radial oscillation of a cloud about 
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its equilibrium position with an exponentially growing amplitude, instead of monotonic 
inward motion. In our case, since clouds can not cool down below a certain temperature 
(equations 7,8), overstability is unavoidable. At each radius a cloud would oscillate about 
its equilibrium radial position with its local buoyant oscillation (Brunt-Vaisala) frequency, 
uj bv = ^(§^7 — f^jff ), where fl c is orbital frequency. We suggest that such radial 
oscillations would provide a natural mechanism to excite acoustic and gravity waves. 
Heating of cooling flow regions by acoustic waves has been considered before (Pringle 1989) 
and will not be considered further here. We will instead focus on the internal gravity waves. 

In cooling flow regions, T and p each may be assumed to bear a power-law form, 
T oc r 7T and p oc r 7p , and if the dark matter density is assumed to be pr>M r lDM , then the 
Brunt-Vaisala frequency uj bv oc t 1dm {^t — §7p) with a positive front coefficient. For normal 
X-ray clusters, 7dm < 0, 7t > 0, 7 P < 0, thus, ojbv is positive and increases monotonically 
towards the center of the cluster. Gravity waves with frequency u > ujbv are evanescent. 
Therefore, in our case where oscillatory motions of cooling clouds are responsible for 
exciting the internal gravity waves, the fact that the buoyant oscillation frequency increases 
towards the center of the cluster suggests that any internal gravity waves generated by 
radially oscillating clouds at a certain radius will be trapped to the region interior to that 
radius. These trapped internal gravity waves will eventually dissipate to heat the gas. 
When considering wave excitement by orbiting galaxies (Miller 1986), Balbus & Soker 
(1990) and Lufkin, Balbus, & Hawley (1995) pointed out earlier that the central regions 
of cooling flow clusters are capable of trapping finite-amplitude resonant internal gravity 
waves, when the orbiting frequency of galaxies falls below the Brunt-Vaisala frequency at a 
certain radius. However, Lufkin et al. (1995) concluded that internal gravity waves excited 
by orbiting galaxies may not be able to provide enough energy to balance the cooling in the 
cooling flow region. 

The energy source for powering the oscillations of clouds hence the excited internal 
gravity waves in our case is gravitational potential energy of cooling gas. Tapping into the 
gravitational energy to heat up cooling flows is not a new idea and its desirable property 
of being able to provide overall balanced heating has been recognized (e.g., Fabian 2003). 
Here, gravitational energy released by a cloud originating at some radius is telescopically 
deposited within some smaller radius at which it finds its quasi-equilibrium position in the 
comoving flow. Heating is hence ultimately due to gravitational energy released as a direct 
consequence of cooling and inward drift of clouds as well as the slower and subsonic inward 
motion of the general flow. Since heating is in essence "reactive" to cooling, cooling should 
be balanced by heating over time, although some spatial and temporal fluctuations may 
be expected. An energetic advantage of heating by the trapped gravity waves is that it is 
economical without significant waste because of trapping. 
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Another interesting property of this telescopic heating by gravity waves is that all 
excited waves at radius larger than r would contribute to heating of gas at r. A somewhat 
more quantitative and relevant illustration of heating by internal gravity waves may be 
made as follows. Assuming mass dropout rate per unit radius is constant, as observed, 
i.e., dM/dr = C, in the cooling region r = — R, assuming that the temperature scales 
with radius as T(r) = T (r/i?) 7T and temperature traces local gravitational potential, 
kT(r)/m p = GM(< r)/2r, and gas density goes as p(r) = po(r/R) lp , assuming that 
interval gravity wave energy is deposited uniformly in radius interior to its trapping radius, 
assuming that each cloud's gravitational potential energy going from its initial position (r-j) 
to the final position (77), = f* f g(r)dr, is released interior to its equilibrium position, 
assuming that CR pressure evolves adiabatically from r\ to 77, and assuming that all gas 
clouds that move inward cool to the allowed floor temperature, then a simple calculation 
yields the gravity wave energy deposition rate interior to a radius r to be equal to 



E(< r) 



R 



-C 

x m 



dx + 



/ 

Jo 



" C^tdx 



2M(< r)kT(r) 
rriplT 



1 + P 



37t/(7p~37t) 




1 + 

37t/(7p~37t) 



1 + 7t/ \l + a~\ 



3/(7p~37t) 



- 1] 



where r\ is the initial position of clouds that find their equilibrium position at r. The 
normally implied heating rate, assumed to be equal to the mass dropout rate (assuming 
heating balances cooling), would be E(< r)\ local , implied = ^ ^^dx = 3 ^ffffi - For 
a = 0.3 and (3 = 0.3, which give T p = 0.33T a , and using 7 T = 0.75 and 7 P = —0.75, 
consistent with some observed X-ray cluster cooling flows (e.g., Kaastra et al. 2004), we 
evaluate the ratio E{< r)/E{< r)\i oca i timp i ied = 48 for R/r = 10, which is comparable to 
the energy dropout rate at R. Therefore, it is seen that the telescopic heating results in a 
progressively stronger heating rate at small radii than would be implied by its local observed 
mass dropout rate. This unique property may help explain the large inferred/required 
heating rate of H a emission-line nebulae (Heckman et al. 1989). But we caution this simple 
calculation serves only an illustrative purpose. 

A desirable property of this heating mechanism by gravity waves is self-regulation, as 
insightfully pointed out by Balbus & Soker (1990). If, say, the central region is overheated 
such that the central temperature becomes very high and possibly the central temperature 
inversion is removed, that would then cause the Brunt- Vaisala frequency to drop towards 
that central region. Consequently, inward propagating gravity waves generated outside 
would get reflected at the surface enclosing the overheated region which marks the peak 
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of the Brunt-Vaisala frequency. Thus, the large amount of heating provided by waves 
generated at large r which caused the overheating in the central region is no longer able 
to deliver energy there; the heating rate of the overheated region would dramatically 
decrease to be unable to balance the cooling and the temperature would drop, until a rising 
Brunt-Vaisala frequency towards the center is achieved. 

Finally, it might be expected that trapped quasi-spherical gravity waves perhaps 
eventually lead to quasi-spherical dissipation patterns in the cooling regions with the 
patterns likely being more visible in the central regions due to more intense cooling/heating, 
as quantified above. This feature might have already been seen in the Perseus cluster from 
recent Chandra observations by Fabian et al. (2003). 



3. Discussion 

The critical assumption made is that CR are frozen to thermal gas on the cooling 
time scale. The diffusion time scale of CR is highly uncertain in the absence of reliable 
knowledge of the topology of magnetic field and turbulence. Nonetheless, it may still be 
useful to have some estimates, but bearing in mind that these are just estimates with 
large uncertainties. For pitch-angle scattering diffusion for Alfven turbulence one finds 
the diffusion coefficient along the magnetic field lines (e.g., by integrating equation 15a of 
Dermer, Miller & Li 1996) to be k„ = f ((? _ 1)(2 \ )(4 _ g) ^LC 1 ^^) 2 "^ 2 ^, where v is CR 
velocity, k m i n the minimum wavenumber of turbulence, the non-relativistic Larmor radius 
of CR, £ the normalized energy density in either shear Alfven or fast-mode waves, p the 
dimensionless CR momentum v/c^Jl — v 2 /c 2 and — q the power spectrum index of Alfven 
turbulence. The actual diffusion coefficient depends quite sensitively on the power spectrum 
index of Alfven turbulence, q. For q = (5/3,3/2,4/3) we find, respectively, ku = 1.3 x 
10 29 (f)(T^)- 1/3 (fe) 2/3 CV /3 cm 2 s- 1 , K ||=4.7xl0 27 (H)( T |,)^/2 ( A f ^ ) i/2 C y/2 cm 2 s -i ; 

and = 7.2 x 10 25 (f )(y^)~ 2/3 ( :? ^) 1/3 £ i ~ V /3 cm 2 s _1 , where \ max is the maximum 
wavelength of turbulence. CR diffusion across magnetic field lines is generally thought to 
be much smaller and neglected here. Thus, during a cooling time of order 10 8 yr, CR may 
diffuse by a distance of order lOkpc, 3kpc and 0.3kpc, for the three cases. We see that the 
assumption of CR being frozen with thermal gas may be valid for parcels larger than from 
~ 0.3kpc to lOkpc, but it is entirely possible that the diffusion scale could be significantly 
smaller than that range. This issue clearly deserves further attention, most likely to be 
resolved only with better simulations based on first principles. 

The combined thermal-dynamical process is influenced by four time scales: cooling 
time, the oscillation period (the inverse of the Brunt-Vaisala frequency), CR diffusion time 
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and cloud sound crossing time. The implicit requirement for isobaric evolution is that 
cooling time is longer than sound crossing time of a gas parcel. For a parcel of size of lkpc, 
the sound crossing time for a sound speed of 300km/s is 3 x 10 6 yr, whereas the cooling time 
is > 10 8 — 10 9 yr. Thus, for parcels of interest here, it appears that isobaric condition is 
met. The oscillation period is in the range of 10 7 — 10 8 yrs in the r = 10 — lOOkpc range. 
Our derivation assumes conservatively that the parcel of gas is able to cool rapidly enough 
to reach its possible floor temperature before it reaches its radial equilibrium position. 
For a cloud initially only slightly overdense the ratio of the final equilibrium radius to 
its initial radius is 0.39 for a = 0.3 and (3 = 0.3, if cooling is rapid enough. The time of 
the travel to the final equilibrium position should be comparable to or longer than the 
oscillation period or the cooling time, whichever is longer. The possibility that cooling time 
is longer than the oscillation time indicates that the cloud can only fall on the cooling 
time, partially supported by buoyant force on its way in, and that the inward motion may 
not be monotonic. It may be argued that ultimately CR would diffuse out, since diffusion 
time, however long, is likely to be finite. Weak diffusion of CR may in fact be necessary 
to ultimately enable gas to move inward. More detailed calculations possibly require joint 
considerations of all the processes involved most likely in a simulation setting. 

Invariably, CR diffusion time must have a certain distribution, depending probably 
primarily on the topology of local magnetic field within a fluid parcel. For example, a fluid 
parcel may contain a set of relatively open magnetic field lines and its CR would be more 
difficult to trap. Such a fluid parcel would then continue to fall inward further than others 
and would cool to a significantly lower temperature, in a fashion customary in the usual 
thermal instability picture. It is instructive to consider a somewhat more general form 
of cooling-heating function ( = ApT^ — rj, and for illustrative purpose we assume equal 
mass heating (i.e., rj is constant; using equal volume heating gives qualitatively similar 
conclusion), for which the thermal instability criterion becomes 

Clearly, for negative 0, equation (9) is always satisfied hence the gas always thermally 
unstable. One way to understand this is that, with negative 0, even isochoric cooling is 
thermally unstable. Therefore, gas with temperature less than about 0.1 — l.OkeV, depending 
on the metallicity of the gas (e.g., Sutherland & Dopita 1993), where Bremsstrahlung 
cooling no longer dominates and negative, is thermally unstable, regardless of the CR 
pressure. Thus, a dense fluid parcel with "leaky" CR may be thermally unstable on its 
way inward. Once its temperature has dropped below about 0.1 — l.OkeV, the runaway 
thermal instability sets in and its temperature may continue to drop rapidly until it reaches 
about 10 4 K, where becomes highly positive and the gas would again become thermally 
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stable. These rapidly cooling clouds may end up near the center of a cooling flow and 
may be identified with emission-line nebulae (Heckman et al. 1989). Since gas parcels with 
higher metallicity would transit to runaway thermal instability at a higher temperature 
in the 0.1 — l.OkeV range, one signature might be that most emission-line nebulae may 
be biased to have higher than average metallicity. Under the same arguments, metals 
may be "transported" in this fashion to the centers of cooling regions, possibly producing 
patches of relatively metal rich gas and unevenness in the metallicity distribution near the 
centers of cooling flow regions as well as creating an overall negative metallicity gradient. 
The expected CR pressure within cold clouds may easily maintain a sound speed of 
100 — 200km/s to explain the observed line width of the nebulae. 

We have neglected thermal conduction in our treatment, which may be justified given 
the requirement of a significant magnetic field to confine CR. The actual distribution of 
gas mass in the range ~ 0.3T a to T a at each radius, where thermal instability plays, would 
depend on many factors including the density fluctuation spectrum and heating function. 
We will not attempt to make further calculation of it here. 

We suggest that the floor temperature at each radius derived in §2.2 is not expected 
to be visibly correlated with other cluster properties, since it is solely determined by shock 
properties, which are not clearly related to any other quantity. It is important to remember 
that the floor temperature should also exist in non-cooling clusters, such as Coma, as 
observed (Kaastra et al. 2004), since a significant CR pressure should exist in all clusters, if 
produced by cluster shocks. 

4. Conclusions 

We put forth a new model for cooling flows in X-ray clusters of galaxies. The model is 
based on the adoption of the view that a significant fraction of the ICM pressure is in cosmic 
rays, as simulations have shown (Miniati et al. 2001; Ryu et al. 2003) and observations have 
indicated (e.g., Blasi 1999). We show that a significant CR pressure changes qualitatively 
the nature of the thermal evolution of the cooling gas. Specifically, the conventional isobaric 
thermal instability condition is only met for a finite range of gas temperature, while a fluid 
parcel becomes thermally stable when its thermal pressure drops below a certain fraction of 
its CR pressure. As a result, the lowest possible temperature at any radius for bulk of the 
gas is close to one third of the ambient temperature. This explains naturally the observed 
puzzling fact (Peterson et al. 2003; Kaastra et al. 2004), without the need to make any fine 
tuning. 
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Additionally, we suggest that dissipation of internal gravity waves, excited by radial 
oscillatory motions of inward drifting cooling clouds, may be responsible for the overall 
heating to balance cooling. The gravitational energy of inward moving clouds is converted 
to the form of kinetic energy of the clouds, which then oscillate and transfer their kinetic 
energy to gravity wave energy, which in turn is trapped interior to the location of each 
wave generator. Therefore, the ultimate energy source for powering the cooling X-ray 
luminosity and heating up cooling gas is gravitational. The proposed heating process has 
several desirable properties. First, it is relatively widespread spatially. Second, heating 
is spatially distributed and globally energetically matched with cooling. Finally, heating 
is strongly peaked at the center, apparently capable of providing adequate heating for 
the observed emission-line nebulae. Possible signatures of such heating include occasional 
quasi-spherical density enhancements, mostly likely seen in the inner cooling flow regions 
due to an expected rapid drop-off of enhancement with increasing radius. 

This research is supported in part by grants AST-0206299 and NAG5- 13381. I thank 
Hui Li and Jim Stone for useful discussion. 
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